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Abstract 

The shot noise of current through a metallic double quantum dot structure exhibit- 
ing negative differential conductance is studied. We can exactly solve the master 
equation and derive an analytical expression of the spectral density of current fluc- 
tuations as a function of frequency in the first Coulomb staircase region. For a 
large range of bias voltage the noise is calculated by Monte-Carlo simulation. Wc 
show that the noise is always sub-Poissonian though it is considerably enhanced in 
the negative differential conductance regime. 
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The shot noise (SN) is a striking consequence of charge quantization and its study has 
become an emerging topic in physics of nano-devices because measurement of the SN can 
reveal more information on transport properties which are not available through the con- 
ductance measurement alone [1,2]. In the case of uncorrelated current the SN has the full 
(or Poissonian) value 2eJ, where e is the elementary charge and / is the average current. 
Deviations (either suppression or enhancement) from this value are due to correlations in 
the motion of charge carriers. The measure of these deviations is the Fano factor F which 
is defined as the ratio of the actual noise spectral density to the full SN-value. There are 
mainly two kinds of correlation: the Pauli exclusion principle and Coulomb interaction. 
While the former correlation always causes a suppression of SN, the latter may suppress 
or enhance the noise depending on the conduction regime. The non-Poissonian noise has 
been most extensively studied in double barrier resonant tunneling diodes, where the SN 
is partially suppressed in the positive differential conductance (PDC) regime and becomes 
super- Poissonian in the negative differential conductance (NDC) regime [3]. The super- 
Poissonian noise accompanied by an NDC has been also predicted and observed in quantum 
dot devices [4,5]. For Coulomb blockade metallic structures when the Coulomb interaction 
is extremely important and manifested explicitly in the so-called charging energy, the SN 
has been studied in a number of works focussing on the single-electron transistor (SET) in 
PDC regime [6-8] (and refs. in [1]). It was shown that the noise is generally suppressed with 
a magnitude depending on device parameters and on the range of applied voltage. Recently, 
we were able to derive an analytical expression for current- volt age (I-V) characteristics and 
a condition for observing an NDC in a Metallic Double Quantum Dot Structure (MDQDS) 
in the first Coulomb staircase region [9]. The NDC has been analyzed in detail with respect 
to device parameters as well as to the temperature and the off-set charge. Then, it is natural 
to raise a question about the SN-behavior in this device in an NDC regime. We will show 
in the present letter that the noise though considerably enhanced in NDC regime seems to 
be always sub-Poissonian. This has been done by solving exactly the master equation in the 
first Coulomb staircase region and by Monte-Carlo (MC) simulation in a large range of bias 
voltage. 

The equivalent circuit diagram of the structure studied is drawn in the inset of Fig. la. 
Within the framework of the Orthodox theory the state \i > of the system is described by 
the probability p(i,t) to have rij (and rrii) excess electrons in the dot D 1 (and D 2 ). This 
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probability obeys the master equation which can be written in the matrix form as dp (t) jdt — 
Mp (t), where p (t) is the matrix of elements p(i, t) and M is the evolution matrix of elements 
Mij = T(j <— i) — SijJ2kF(k ±— i), with F(j <— i) being the net transition rate from \i > 
to \j > (see [6,7] and formula (1) in ref. [9] for the model under study). The statistically 
averaged current across the junction v {y — I, m or r) is (I u (t)) = e E« [r* (i) ~ (i)] p(i, t), 
where T^(i) is the tunneling rate through the junction v to the right (+) /left (-) at the state 
\i >. For a stationary state the current is t-independent and the total external current is 
/ = (I u ) for any of v. However, if the state is not stationary, the charge accumulated in the 
dots is time dependent and (I(t)) is a weighted average of {h(t)) as (I (t)) = Y. v 9u {Iu {t)), 
where g t = C r C m jYi c , g m = CiC r /E c , g r = CiC m /T, c and S c = CiC m + C m C r +CiC r . We will 
be interested in the SN for J, which is the external current and experimentally measurable 
[10], as well as for partial currents l v . 

Korotkov [7] suggested in detail the procedure of calculating the noise spectrum of a 
correlated tunneling current in a SET. Extending this procedure for the MDQDS of interest 
we have 

S vv (u) = 2A U + 4e 2 E [r+ (i) - V- (i)] B tj [r+ (j \v~ ) p st (j — F~ (j \u+ ) p st (j \v + )\, 
S u M = 2 E 9lA v + 4e 2 E E 9^ [r+ (i) - T~ B l3 (1) 

V U^L ij 

x "rj (j \ir ) Pst {j \n~ ) - r- (j ) p st {j \fi + ) . 

Here S uu and Su are the spectral densities of current fluctuations (SDCF) for currents I v 
and /, respectively; A v = e(I+ + I~) with l^- = eJ2iPst(i)F^; the conditional probability 
p{% <— j\r) for having state \i > at the time t = r > under the condition that the state was 
\j > at an earlier time t — obeys the same master equation as for the probability p(i, i); the 
stationary probability p s t(i) is defined asp(i <— j\r — > oo) = p s t(i)Sij] B = Re{{iuil— M)^ 1 }; 
(jl^) is the state obtained from the state \j >= (rij,mj) by transferring an electron across 
the ^-junction to the right (+) /left (-). 

For the structure under study, using the well-known expression of the tunneling rate 
across a junction (see, eqs.(4,5) in [9]), in principle, one can calculate the SDCFs (1). In 
practice, however, one can not solve the master equation exactly with all possible states. 
Recently [9], we have shown that at zero temperature and in the first Coulomb staircase 
region, V sl < V < V s2 , where V s2 = e/2C r and V sl is the maximum from e/2Ci and 
e\Ci — C rn \/2C r (Ci + C m ) (assuming C\ > C r ), the master equation can be exactly solved 



and therefore the I-V characteristics can be derived for the range of parameters as 

C r < Ci < 3C r and C m < C r (3C r - C7,)/(C, - C r ). (2) 

Under this condition all probabilities p s t(i) are equal to zero except those for three states 
|1 >= (0,0), |2 >= (1,0) and |3 >= (1,-1): p st (l) = 6c/£ r ; p st {2) = ca/£ r and p st {3) = 
ab/T,r, where 



a = T (2 <- 


- 1) = 


[Q (C m + C r ) /eXcRt] (V - e/2Q) , 




b = T{3^ 


-2) = 


[C r (C m + Cj) /e£ c it> r ] (V - e (C, - C m ) /2C r (C, + C m )) , 


(3) 


C = r(i^ 


-3) = 


[aC r /eE CJ R m ] (e/2C, + e/2C r - V) , 





S r = ab + 6c + ca. With (3) the SDCFs (1) become 

S u = 2el (1 + 2aB 12 ) , S mm = 2el (1 + 2a£ 31 ) , S rr = 2el (1 + 2aB 23 ) , 
5// = 2el Egl + 4el[ gi a (g m B u + g t B 12 + ^ r Si 3 ) + {g rn B 21 + ^ + # rJ B 23 ) (4) 
+ 3mC (^ m -B 3 i + 9iB 32 + # rJ B 33 ) . 

The matrix {iul — M) has then a simple form giving the matrix B with elements 

Bij = Re{D i:j [iu;(ab + bc + ca- to 2 ) - uj 2 (a + b + c)] -1 }, (5) 

where D u = be — uj 2 + iu{b + c); Di 2 = 6c; D 13 = be + iuc; D 2X = ac + iua; D 22 = 
ac — oj 2 + ?cj(a + c); _D 23 = ac; D 3 i = ao; _D 32 = ab + iub; D 33 = a6 — cj 2 + icj(a + 6). The 
expression (4) (with (3) and (5)) is our main analytical result obtained for the first Coulomb 
staircase region in the condition (2). Substituting (3) and (5) into (4) we see that for the 
model under study all S vu {uj) (v = l,m,r) are identical 

S V u M _ 2afec (a + b + c) 

2e/ ~ (ab + bc + ca- uo 2 ) 2 + uo 2 (a + b + c) 2 ' 

In the limit of zero frequency the noises SW(0) and S7/(0) are coincident with a single Fano 
factor 

F = [I /a 2 + 1/6 2 + l/c 2 ][l/a + 1/6 + l/c]" 2 . (7) 
In the opposite limit of large frequency F vv — > 1 (the Poissonian value), whereas Fu — > 

For given values of device parameters (capacitances and tunneling resistances) satisfying 
the condition (2) it is easy to calculate the SDCFs (4) in a large range of frequency. To this 
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end we use the zero temperature tunneling rate across a junction T = Q(—AF)\AF\/e 2 R t , 
where R t is the junction tunneling resistance, AF is the change in the free energy F of 
the system after the tunneling event has occurred. For the model under study F(i) = 
(eni - CiV/2f/2Cl + (errii + C r V /2f /2C* r + {en { - Q\//2)(em i + C r V/2)/2C* m + eV(m - 
n r )/2 - {d + C r )V 2 /8, where C* = E C /(C, + C m ),C* r = E C /(C, + C m ),C* m = Y, c /C m ; 
ni(n r ) is the number of electrons that have entered the structure from the left (right) [9]. In 
calculations as well as in MC-simulations below the elementary charge e, the capacitance C r 
and the tunneling resistance R r are chosen as the basic units. The voltage, the current and 
the frequency are then measured in the units of e/C r , e/C r R r and (C r i? r ) _1 , respectively. 
Fig. 16 shows the normalized SDCFs S uv {u) /2el (two upper lines) and Su(u>)/2el (two lower 
lines) calculated at the same bias voltage V = 0.5 for two cases of I-V characteristics with an 
PDC (dashed line) and an NDC (solid line) as shown correspondingly in Fig. la. This result 
should be in comparison with that for the SET shown in Fig.3 of [7]. For the more interesting 
case of NDC the Fano factor is plotted versus the bias voltage in the inset of this figure. 
Clearly, the noise is considerably greater in the NDC regime, the factor F is however still 
limited by the Poissonian value. An enhancement of the SN in an NDC regime is observed 
in various structures due to a strong fluctuation of current. The present result of F < 1 
is different from that observed in resonant tunneling diode devices where the electrostatic 
potential fluctuation of the band bottom leads to a super- Poissonian noise [1,4]. It should be 
emphasized that the origin of the super-Poissonian noise is related to the nature of potential 
fluctuations producing a positive feedback of charge which is absent in the MDQDS under 
study. On the other hand, the sub-Poissonian noise accompanied by an NDC has been 
observed in GaAs — AlAs — GaAs tunneling structure with embedded self-assembled InAs 
quantum dots in the single-electron tunneling regime [11] and also suggested in strongly 
correlated double quantum dot systems in the Kondo regime [12]. 

For a large range of bias voltage we simulate the noise using basically the MC-program 
[9] which was shown to give I-V curves in good agreement with analytical calculations (see 
Fig. 2 in [9]). However, it is impossible to simulate the spectral density in the limit of zero 
frequency and we have to consider the low-frequency limit uo c = 10 -3 (indicated by the 
arrow in Fig. 16). For all frequencies u > oo c the simulation noises are practically coincident 
with corresponding analytical curves in Fig. 16. As an additional test, our noise program 
has well reproduced the qualitative behavior of experimental data in Fig.l of [2] for a single 
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dot structure. Fig. 2 presents the voltage dependence of normalized simulation SDCFs for 
sample with parameters given in the figure. The calculation has been performed at uj = uj c 
in a range of V where the I-V curve (dashed line) exhibits several peaks and valleys. Clearly, 
two solid curves describing S UU (V) and Sjj(V) have practically the same form, though at 
the chosen finite frequency there is still a considerable separation between them. At each 
bias voltage the Fano factor F is somewhere between these normalized noises and we can 
guess that the F(V)-dependence should have the same saw-tooth behavior like solid curves 
in Fig. 2. The most interesting feature observed in this figure is the modulation of the 
noise amplitude as a function of V with peaks at the points where the NDC reaches the 
highest magnitude and moreover the noise seems to be always sub-Poissonian even in NDC 
regimes. Such a suppression of noise is due to the strong Coulomb interaction as generally 
suggested in ref.[8]. We like to mention that a similar voltage-dependent behavior of the 
Fano factor has been experimentally observed in [11] for a single quantum dot structure. The 
gradual decrease of noise peaks at NDC regimes as the voltage increases is closely related 
to the corresponding decrease of the peak-to- valley ratio of current as can be seen in Fig.2. 
Note that qualitatively all the simulation results discussed are not particular for the chosen 
frequency iv c . Thus, both the analytical results of eqs.(6) and (7) and the MC-simulation 
data suggest that the SN in the MDQDS under study is always sub-Poissonian though 
it is considerably enhanced in NDC regimes. This is the main conclusion of the present 
work which, as shown by additional simulation data (should be published elsewhere), is well 
preserved even in the case gates are included. 

One of authors (P.D.) acknowledges the CNRS for financial support under PICS No. 404. 
The work in Hanoi was supported in part by the Natural Science Council of Vietnam and 
by the VAST research grant on simulation of nano-devices. 
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FIGURE CAPTIONS 



FIG.l The I-V characteristics (a) and corresponding normalized SDCFs calculated at V 
= 0.5 (b) for Ri = 1.1 (PDC - dashed lines) and Ri = 0.2 (NDC - solid lines), everywhere 
Ci = 1.5, C m = 2.0, R m = 2.0. Insets: in (a): circuit diagram of the model, in (b): the Fano 
factor as a function of the bias for the NDC case. 

FIG. 2 The voltage dependence of normalized SDCFs: solid curves (upper for S uu and 
lower for S n ) at the frequency u c = 10~ 3 . The corresponding I-V characteristics is shown 
by the dashed curve. Simulation parameters: C\ = 1.0, C m = 1.0, Ri = 1.0, R m = 10.0, zero 
temperature. 
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